The most generalized analytical approximation to the solution of single-mode 
spin-boson model without rotating-wave approximation 

T. Liui'2, K. L. Wang3, and M. Feng^* 
State Key Laboratory of Magnetic Resonance and Atomic and Molecular Physics, 
Wuhan Institute of Physics and Mathematics, Chinese Academy of Sciences, Wuhan, 430071, China 
^ The School of Science, Southwest University of Science and Technology, Mianyang 621010, China 
The Department of Modern Physics, University of Science and Technology of China, Hefei 230026, China 

The single-mode spin-boson model (SMSBM) has extensive application in different subfields of 
physics. In the absence of rotating-wave approximation (RWA), we try to solve SMSBM analytically. 
We argue that the analytical expression obtained is the most exact approximation to the solution 
of the system under the assumption of Abel-Ruffini theorem, which works for a wide range of the 
parameters such as coupling strength and detuning and would be practical for currently available 
experiments. 
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The single- mode spin-boson model (SMSBM), describing a two-level system experiencing a single- mode boson 
field, is an important prototype in diverse phenomena in almost every subfield of physics. From early day's studies 
of Holstcin model [1] in condensed matter physics and Jayncs-Cummings model [2] in quantum optics to recent 
investigation of quantum information processing [3, 4], SMSBM has been playing crucial roles. 

Rotating-wave approximation (RWA) has been usually employed in treating SMSBM to simplify the solution under 
the condition of near resonance and weak coupling within the characteristic time of the system. However, RWA is not 
working well any more recently due to experimental availability of strong coupling in atomic system [5] , semiconducting 
system [6] and superconducting system [7], which imphes the necessity of solution to SMSBM in the absence of RWA. 

Without RWA, however, it is hard to have an analytical solution to SMSBM [8] because the counter-rotating terms 
make the computational subspacc unclosed. As a result, no matter what methods were taken [9, 10, 11, 12, 13, 14, 
15], numerically or scmi-analytically, the solutions were made based on the truncated subspacc under some special 
conditions. Alternatively, as the coherent state consists of infinite numbers of Fock states, the computational subspace 
in the absence of RWA, which is unclosed in the basis of Fock states, could be nearly closed in the coherent-state 
representation for a wide parameter range. Therefore, employment of coherent states would perhaps enable us to have 
analytical approximations very close to the exact solution of the problem. 

There have been some peculiar characters discovered in SMSBM without RWA, such as Bloch-Sicgcrt shift [16], i.e., 
a shift with respect to the true resonance frequency due to counter-rotating terms, and quantum chaos in cavity QED 
by differently polarized lights [12]. It has also been shown that the discrepancy of the SMSBM in the presence of RWA 
with respect to the absence of RWA is reflected by some phase dependent effects [13]. Most of those discrepancies are 
meaningful only theoretically, whereas Bloch-Siegert shift is observable experimentally. For example, in ion traps, the 
strong laser radiation on the trapped ultracold ions could lead to some level deviation regarding Bloch-Siegert shift 
[17, 18, 19]. As the laser is usually treated as external classical light, counter-rotating terms could show observable 
effects in blue detuning case. Recent study has demonstrated the possibility of fast logic gating with strong coupling 
between laser and trappcd-ion qubits, in which the usual RWA treatment could not work well [20], and the strong 
nonlinearity would yield considerable complexity in the time evolution [21, 22, 23]. As a result, it is of importance to 
have a more strict study of the SMSBM without RWA. 

We have noticed a very recent publication for analytically solving SMSBM by generalizing RWA [24]. In the 
present paper, we will try to explore a more efficient approach to a more generalized analytical result, compared to 
[24]. The key point is that, by employing coherent-state representation, we try to diagonilize the matrix regarding 
the computational subspace. Based on the Abel-Ruffini theorem that no general solution in radicals is possible to 
polynomial equations of degree five or higher [25] , we present analytically the expressions of the energy levels of the 
system under third-order approximation of the exact wavcfunction, which is relevant to a polynomial equation of 
degree four. We argue that this should be the most generalized approximation to the solution of the SMSBM. To 
further show the validity of our solution, we will give evidences both analytically and numerically in comparison with 
the results in previous publications. Moreover, we will show potential application of our result. 
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Consider following Hamiltonian in units of ?i = 1 [24] , 

f u;oa+a + \{a+ + a) \ , . 
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where is the energy splitting of the spin and wq is the frequency of the boson field. A denotes the coupling between 
the spin and the boson field, and a (a+) is the annihilation (creation) operator of the boson field. Eq. (1) is taken 
from [24] in order for us to make comparison between our solution and in [24]. In fact, with respect to the standard 
quantum optical notation, Eq. (1) has been taken a unitary rotation on the two-level system. However, as there is 
no change of physical essence with that unitary transformation, we will work on Eq. (1) in most of this paper. 

For convenience of treatment, we may set g = A/wq a-nd employ displacement operator D{g) = exp [17(0^ — a)] 
acting on and a, which yields A = D{gy aD{g) ~ a + g, = D{g)'^ D{g) ~ + g, B ^ D{—g)''aD{—g) = a — g, 
and = D{—g)'^a^D{—g) = at — g. Consequently, Eq. (1) could be rewritten as 

ojo{A+A^g'^) n/2 
n/2 LJa{B+B- 

which is formally solvable, and we assume following trial solution to Eq. (2), 

N 
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where ]e) — and \g) ~ ^j*), c„ and (i„ arc coefficients determined later, and is a large integer relevant to 
the size of the truncated subspaccs. |?t.)a(_b) is a coherent state regarding the operator A{B), which defines as 
\n)A = -^(flt + g)"lO)^ and \n)B ~ "^("^^ ^ .9)"|0)-B with |0)^(b) the coherent state in the subspacc regarding the 
operator A{B) [18]. 

Putting Eq. (3) into the Schrodinger equation of Eq. (2) could yield 

uja{m - g'^)cm + -^^^{-iTD^^d^ = Ecm, (4) 
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[18]. Eqs. 



where (-!)"£)„„ = (-l)™D„i„ = s(mln)A, and Z3mn = e I]fe=o (m-fc)!(n-fc)!fc! . . 

(4) and (5) present the possibility to have a closed solution to the problem. To analytically solve Eqs. (4) and (5), we 



set dn ~ ±(— l)"c„, which yields tooim 
on following determinant, 
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where e^j = cJo("^ — .9^) + ^m.m ^ = ^ with il^j „ = ±(l/2)r2Dm„. As the superscripts ± are consistent 
for the relevant variables, the cases regarding superscripts '+' and '— ' will be treated independently. In principle, if 
we consider a large enough value of TV, Eq. (6) would lead us to a nearly exact solution to Eqs. (4) and (5). However, 
in terms of Abcl-Ruffini theorem, to have the analytical results to the best, we have to reduce the determinant to 
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FIG. 1: Energy levels with respect to (A/ljo), where we assume (a) wo = and (b) = 0.750. The solid curve and circles, 
respectively, mean the numerical (exact) solution and the solution from Eqs. (8) and (9) under the third-order approximation. 

which leads to a polynomial equation with degree four and should be the most generalized analytical approximation 
of the solution to the SMSBM under our consideration. Straightforward deduction to Eq. (7) yields 
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with m = 0, 1 corresponding to the ground and the first excited states, respectively. Other excited states are 
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To prove the validity of Eqs. (8) and (9) explicitly, we may work along following two aspects: Comparison with 
numerical treatment of Eq. (6) in the case of a big enough value of N, and comparison with other analytical solutions 
by the determinants with smaller subspaces. For the former, we have made numerics on Eq. (6), as shown in Fig. 1, 
by setting iV = 42 with off-diagonal elements fiij {i or j > N) smaller than 10~^. We may consider that numerical 
result is the exact solution to the problem. Fig. 1 shows that our results in Eqs. (8) and (9) under the third-order 
approximation agree with the exact solution very well even in the case that A, f2 and luq are comparable. 
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FIG. 2: Energy levels with respect to [X/uq) by Eq. (10) under the first-order approximation (solid curve) and by the results 
from [24] (the circles), where we assume loq — 0.75S1, (a) and (b) correspond to ground state and excited states, respectively. 



For the latter, we first consider the zero-order approximation of Eq. (6), 
(A^/ciJo) ± The first-order approximation corresponds to 



0, 



0, which yields £"„ 



(10) 



which leads to the ground state energy, Eq = c^oi^ ^ 9^) ^ §(-Do,o + ^i.i) 



and the energies for excited states E\ 
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ujQ- {-lYmDk,k- Dk+i,k+iW' +^"^01^^^^, where k = 0,l,---,iV and Dk,k+i = i'fc+i.fc is used. It 

could be found in Fig. 2 by comparison with the results in [24] that the excited-state energies we obtain are in good 
agreement with those in [24], but the ground state not. In this context, we consider that the results obtained in [24] 
is basically the one under the first-order approximation in our treatment. But the ground state plotted in [24] seems 
to be Eq = —Q,Dqx>/2 — A^/ojo, i-c, the zeroth-order approximation shown above [26]. 



Analogically, we may also obtain the second-order approximation using 



'''•m+l.m 
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for which we omit the lengthy expression of the analytical result, but emphasize that the accuracy of the solution 
depends on how many off-diagonal terms in the determinant of Eq. (6) are involved. In general, the farther the 
off-diagonal elements away from the diagonal line of the matrix in Eq. (6), the less significant the elements play their 
roles in the solution. But in the case that O, u>o and A are comparable, our numerics shows that the situation is 
very complicated, for example, fis.o, f^2.o- ^o.Si ^0.2 becoming comparable to fli^ and ^o.i- This is the reason that 
omission of the elements other than the nearest neighbor to the diagonal terms of the matrix yields the deviation in 
the mediate coupling case in [24] with respect to the exact solution. In contrast, our treatment could present results 
more accurate than under the standard RWA and than in [24]. As demonstrated in Fig. 1, our analytical expression 
fits the numerical results very well in a wide range of parameters. Under the assumption of Abel-Ruffini theorem, 
we argue that the results we present in Eqs. (8) and (9) under the third-order approximation should be the most 
accurate analytical expression for the energy levels of the SMSBM under our consideration. 

Why could we make this ? The key reason is the correlation between c„ and d„ we found, i.e., d„ = ±(— l)"c„, in 
the coherent-state representation. The coherent-state representation helps us to have a close subspace for solution in 
the absence of RWA, and the coefficient correlation significantly simplifies our analytical deduction, which makes it 
available to reach the expression under the third-order approximation. 

Compared with purely numerical treatments, our analytical result could present some physics more clearly. For 
example, for n = in Eq. (3), we have the ground state |) = (l/\/2)(|0)/i|e) — |0)s|.g)), which implies that the groimd 
state of the system always overlaps with the upper level of the spin and would always keep evolving if we involve the 
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FIG. 3: Time evolution of the population in the lower level of the trapped ion by Eqs. (8) and (9) (solid curve) and by the 
first-order approximation (dashed curve), where we assume g = 0.8 and the initial vibrational state in coherent state with 
\a\ — 1.0. (a) Qo/i^ = 1 and (b) ^lo/u =3/4. We consider the solid curve is a nearly exact description of the dynamics of the 
trapped ion. As the generalized RWA treatment in [24] is equal to our treatment under the first-order approximation, we may 
find the deviation of the approximate treatment from the exact description with time. 

counter-rotating terms in our treatment. In contrast, under the framework of RWA and even in generalized RWA 
treatment [24], the ground state of the system is always uncoupled from other states and thereby remains unchanged 
no matter how strong the interaction is. Other potential application could also be found in [18]. 

On the other hand, as our analytical result is very close to the exact solution, wc may employ it to study quantum 
behavior of the SMSBM under arbitrary conditions. For example, wc may accurately calculate the dynamics of the 
system in the regimes of mediate and strong coupling, which is helpful in experimentally exploring the dccoherence 
and operational infidelity regarding qubits in quantum information science. Specifically, for strong coupling case in 
trapped ion system, we may unitarily transform the original Hamiltonian to a Hamiltonian very similar to Eq. (1) 
[18]. 

H' = — 2''^^ ^ i^a^a + g{a^ + 0)0^ + ea^ + 5^, 

where fio is the Rabi frequency regarding laser-ion coupling, v is the trap frequency, g is related to Lamb-Dicke 
parameter, and e is the detuning of the laser with respect to the trapped ion. Ox.z are usual Pauli operators based 
on the two levels of the ion [18]. As an example, we demonstrate in Fig. 3 the time evolution of the population 
in the lower level under the initial condition ^^jv = 1 or 3/4, e = and g — 0.8 with coherent state a = 1.0 in 
vibration. As the results based on Eqs. (8) and (9) are nearly the exact solution, we could study the dynamics 
efhciently from the complicated evolution of the system. The figure also shows the deviation of the results by [24] 
from the exact dynamics. What is more, as g = 0.8 implies the case beyond Lamb-Dicke limit [18], our analytical 
results are also useful for understanding the behavior of the ions outside the Lamb-Dicke regime [27] from a purely 
quantum mechanical viewpoint. 

It seems also possible to apply our method to more complicated situation regarding strong spin-boson coupling. 
For example, extending SMSBM to multi-spin singlc-modc boson interaction reaches Dicke model [28]. It has been 
found that the Dickc model in the absence of RWA owns some unique characters [29] with respect to the case in the 
presence of RWA [30] such as in quantum phase transition, in Berry phase and in entanglement. Moreover, considering 
multi-mode boson field, we could have relevance to another fundamental problem with environment interrupting a 
spin-qubit, as shown in a recent work [31] that Zeno effect is stronger in the absence of RWA than in the presence of 
RWA, and anti-Zeno effect disappears if RWA is removed. With our method by minor modification, we may enable 
some analytical discussion for above relevant problems. 

In conclusion, we have presented some analytical expressions for solving the SMSBM in the absence of RWA. We 
argue that our third-order approximation is the most accurate analytical result, which could effectively replace the 
exact numerical solution for studying the SMSBM under arbitrary condition. As SMSBM has been widely applied to 
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different physical problems [8, 32]. the coherent-state forms of the eigenfunction and the analytical expression of the 
energy levels of the system would be helpful for understanding the interaction and the dynamics in the spin-boson 
model under strong coupling or other extreme conditions. 

We also conjecture that our technique could be extended to multi-spin or multi-mode cases. Compared to numerical 
solutions to these cases, our analytical results would help us to get more physical insight from the complexity. 

The work is supported by NNSFC under Grant No. 10774163 and by the funding from the State Key Lab of 
MRAMP. 
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